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Abstract 



We show that the full symbol of the Dirichlet to Neumann map of the A:-form Laplace's 
equation on a Riemannian manifold (of dimension greater than 2) with boundary de- 
termines the full Taylor series, at the boundary, of the metric. This extends the result of 
Lee and Uhlmann for the case k = 0. The proof avoids the computation of the full sym- 
bol by using the calculus of pseudo-differential operators parametrized by a boundary 
normal coordinate and recursively calculating the principal symbol of the difference of 
boundary operators. 



1 Introduction 



While there many resuhs on the uniqueness of recovery of the coefficients of an el- 
hptic partial differential equation from boundary data in the case of a single partial 
differential equation, there are few results for systems of PDEs. One might expect that 
the complete boundary data for a system might be sufficient to recover multiple coeffi- 
cients and yet the results to date have been in essentially scalar cases. Lee and Uhlmann 
showed that the full Taylor series, at the boundary, of a metric can be obtained from 
the total symbol of the Dirichlet to Neumann map of the scalar Laplace's equation. 
One might expect to obtain at least the same information from the Dirichlet to Neu- 
mann map associated with the Laplacian operator on ^-forms. The 1-form, or vector 
Laplacian on 3-manifolds being the example with the most obvious applications. We 
will show here that the full symbol of the A:-form Laplacian does indeed determine the 
Taylor series of the metric at the boundary, and hence under suitable assumptions an 
analytic metric can be recovered from this data. The method used, in common with f4l, 
avoids the computation of the full symbol of the Dirichlet to Neumann map. Rather, by 
using the calculus of pseudo-differential operators parametrized by a boundary normal 
coordinate, we recursively calculate the principal symbol of the difference of boundary 
operators, checking that it vanishes to a suitable order As well as being an natural ex- 
tension of the scalar case and interesting in its own right, we hope that this paper will 
stimulate the use of the technique in other inverse boundary value problems for elliptic 
systems of equations of interest to applications. 

The context for all will be a smooth compact orientable manifold with boundary 
M, equipped with a Riemannian metric g. We also assume dimM = « > 2. The metric 
tensor induces a volume form /j £ Q."{M) and Hodge star isomorphism * : D!^{M) — * 
n"^'^(M) is defined by the property 

*(BA(B = g(a),(X))^, (LI) 

where the action of the metric is extended naturally to act on ^-forms. We can consider 
the Hodge star on ^-forms as a contraction of the tensor ^"^^gj^u. Here is the 
covariant metric tensor. 

The total symbol of an operator P on functions on R" is 

A classical pseudo-differential operator of order m has a full symbol which is an asymp- 
totic sum of terms Pm-j (-x, ^) which are smooth in ^ ^ and for A- > are homogeneous 
of degree m — j 

p,n-jix,Xt,) ^ X'"'^p,„-j{x,t,). 

The principal symbol is p,„ also denoted by a,„(f). The class of classical pseudo- 
differential operators is denoted by M30J,"(M"). There are more general classes of 
pseudo-differential operators based on more general symbols, but we shall not need 
them here. These classes form a graded algebra under composition 

o : TOO^^I X MDO^^f ^ TOO3+"' . 
To obtain the principal symbol of the composite one just takes the product: 
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however the full symbol of the product is rather more complicated. Operators in 
^100^^°° = n ^W)0"j are called smoothing operators. The full symbol of a pseudo- 

mGR 

differential operator determines the operator modulo smoothing operators. For brief 
introduction to pseudo-differential operators we recommend the notes 1 3 1 and for more 
detail Shubin 1101 . We note that the definition of pseudo-differential operators can be 
extended to smooth manifolds using coordinate charts. Here the principal symbol is in- 
variantly defined as a function on the cotangent bundle while the total symbol depends 
on choice of coordinates. 

Following j4l we will consider pseudo-differential operators on a smooth manifold 
Y depending smoothly on a parameter f . For our purposes we will have Y — dM and 
t the normal distance from the boundary. We say that P e *HDO"'''^(y,M+) if it is a 
family of pseudo-differential operators of order m on Y, varying smoothly up to f = 0, 
and such that 

with Pj a smooth family of operators on Y of order m — J. This definition extends 
naturally to operators on bundles, in our case the bundle of A:-forms being the important 
example. 

The symbol ofPe TOO'" ' (y,M+) is defined to be the vector 

K->(^>));=o 

evaluated at f = 0. This is a vector of functions on the cotangent bundle of Y. For the 
case of an operator on a vector bundle, each of these functions is a field of enodmor- 
phisms on the fibres of the bundle. 

Let M be a 1-form then the Bochner Laplacian (sometimes called the rough Lapla- 

cian) is the operator expressed in coordinates as — Y.g''''^k;ij- The principal symbol in 

U 

this case is gl where / is the identity on 1-forms. The formal adjoint with respect to a 
metric of the exterior derivative on A:-forms is 5 = (— 1)^"'^+"+^^ The Laplacian on 
differential forms (or Hodge Laplacian) is then A = dd + 5d. The principal symbol of 
d is a'/(^)a) = A CO. Let X ^03 denote the contraction of the differential form co with 
respect to the vector field X. We denote by the vector field dual to the one form ^. 
The principal symbol of 5 is then CJg(^) — —i^^ -i . We conclude using that contraction 
is an antiderivation on forms 

oi(^)(o = ^ A CO) + ^ A (^« ^ )C0 = ^)co. 

The connection between the Laplacian and the Bochner Laplacian, as well as an alter- 
native way to calculate the principal symbol of the former, is given by the coordinate 
expression for the Laplacian 




(Am; _ 

^a=ip=l / 

Note that for a flat metric the Laplacian and Bochner Laplacian coincide. A differ- 
ential form M satisfying Laplace's equation Am = is called a harmonic form. On a 
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compact manifold without boundary, this is equivalent to the condition that the form is 
a harmonic field, that is it is both exact, du = 0, and co-exact, 5m = as 

(m,Am) = ||£/m|P + ||5m|P+ / 5mA*m+ / uA^du (1.2) 

JdM JdM 

However on manifolds with boundary there can be harmonic forms which are not har- 
monic fields. 

Closely related systems of elliptic partial differential equations occur in electro- 
magnetics (the vector Helmholtz equation) and in linear elasticity. 

In a linear elastic solid with metric tensor g and with no body forces, the displace- 
ment field M (a vector field) satisfies the equation Div(CL„g) = where m L^g is 
the Lie derivative of the metric which in components is (Lug) = Uij + m^;,- (as usual 
a semi-colon indicates covariant differentiation with respect to following indices) and 
Div is its formal adjoint H'^ij kg''^- (^H sums will be indicated explicitly.) The 

elastic tensor C is a field of automorphisms of the symmetric tensors on each fibre. The 
principal symbol of the elastic operator is C. For an isotropic solid C = A,g Cg) + ^/ 
where / is the identity operator on symmetric tensor fields. The problem considered by 
Q was the recovery of the Lame parameters X and fj for an isotropic solid. They also 
considered a related anisotropic problem for two-dimensional elastic media. 

Nakamura and Uhlmann ijs] derive a factorisation for the anisotropic linear elastic- 
ity operator in boundary normal coordinates (for the flat metric). This aflows them to 
recover the full Taylor series of the 'surface impedance tensor', which is a function of 
C, but not the complete Taylor series of C. For a special case, transversely isotropic 
media, C can be recovered 

In electro-magnetic theory the electric field E and magnetic filed H are naturally 
defined as 1 -forms, as to take measurements of these fields one must integrate over 
curves. The resulting electric and magnetic fluxes, D and B are naturally two forms as 
one must integrate them over surfaces to make a measurement. The material properties 
(for simplicity we consider a non-chiral, linear, insulating material) are the permittivity 
£ and permeability /j, these map one forms to two forms and are the Hodge star oper- 
ators for an associated electric and magnetic Riemannian metric. Assuming all fields 
to be time harmonic with angular frequency co and the electric charge density to be 
constant we have Maxwell's equations 

dB =0, <iD = (1.3) 
dE =-m^H,D = &E (1.4) 
dH =mEE,B=ijH (1.5) 

In a physically artificial situation where /j = e = * (obviously after units have been 
scaled) we notice that E and H salsify the vector Helmholtz equations AE = (ir'E and 
A// = co2//. 

The main result we prove is an extension to A:-forms for A: > 0, of the result of Lee 
and Uhlmann |6| using a similar factorization in boundary normal coordinates. Ini- 
tially, our notion of a Dirichlet-to-Neumann map is non-standard for ky^Q. Employing 

the multi-index notation / = (/i , . . . , i^) we write a A--form as m = J^uidxj where dxj = 

I 

dxi^ A c/x;, A • • • A dxii,. Following 16,| we use a coordinate chart {xi,... ,x„) = {x/ ,x„) 
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where x' — (xi,... ,x„^ i ) is a chart on the boundary, and x„ is the distance to the bound- 
ary. We denote by d„ both a unit vector field normal to the boundary and its associated 
derivation on functions. We extend this to k forms by its actions on the components as 
functions d„u = Y,d„ujdxj. The operator '.uldM ^ {'^nu)\dM^ where Au = 0, is linear 

while somewhat unnatural. We will discuss the relationship to more natural Neumann 
data for Laplace's equation in Section 2. Our main result is 

Theorem 1.1. Let M be a smooth compact orientable Riemannian manifold with bound- 
ary, with dim(M) ~n>2 and metric g, and let k be an integer < k <n. 

(i) The Dirichlet-to-Neumannmap A.g ; m|3m ^ {d„u)\^!^^ forthe k-fonn Laplace's equa- 
tion Am = is a classical pseudo-differential operator of order one. 

( ii) The Taylor series, at the boundary, of the metric in boundary normal coordinates 
is uniquely determined by the full symbol of Kg. For < k < n only one diagonal 
component of the full symbol is needed corresponding to dxj — dxi^ Adxj^ A • • • Adxjj^ 
but ft? r k = (n + 1 )/2 the multi-index I — (i'l , . . . , 4) must exclude n and for k — [n — 
1 ) /2, / must include n. 

Where our work differs from Lee and Uhlmann's is in its use of families of op- 
erators parameterized by the normal distance. The case A: = « is clearly equivalent to 
= so we need only consider the case < k < n. Lee and Uhlmann showed that 
the full Taylor series of a real analytic metric on a real analytic manifold, where the 
relative homotopy group of the boundary K{M,dM) is trivial, determines the metric, 
provided the manifold is strongly convex or the metric can be extended analytically to 
a larger manifold without boundary. Recent work of Lassas and Uhlmann |5 j removes 
these geometric and topological hypotheses showing that an analytic metric is deter- 
mined throughout any connected analytic manifold by the Dirichlet-to-Neumann map. 
It remains to be seen if the same techniques can be applied to the general A:-form case. 

Only a small modification of the argument is needed to prove a version of Theo- 
rem lLll for the equivalent A;-form Helmholtz problem at fixed frequency. 

2 Boundary Conditions 

In a neighbourhood of the boundary where our boundary normal coordinates are de- 
fined, we can distinguish between tangential ^-forms which have no dx^ in their co- 
ordinate expression, and normal forms which have a common factor dx^. Clearly the 
space of A;-forms on this neighbourhood is the direct sum of the spaces of tangential and 
normal forms. The projection on to the tangential component is 7it(co) = 3„ {dx,, A Co) 
and on to the normal component Tin (to) = dx„ A (3„ ^co). From these expressions it is 
clear that *Ka = K^* and that the Hodge star of a tangential form is normal and vice- 
versa. 

Let i : dM M be the inclusion of the boundary. The tangential component 
of ^-form u S n^(M) at the boundary is then the pull-pack to the boundary i*u G 
Q.^{dM). The normal part of u at the boundary, can be determined uniquely from 
i* *u = *a(3n -'M)|aM G Q."^^{dM) where *^ is the induced Hodge-star on the bound- 
ary. Dirichlet data for harmonic fc-forms consists of the both the tangential component 
and the normal component of the form at the boundary (T|. Note a possible source of 
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confusion here. When considering forms on manifolds with a boundary i*u is often 
thought of as the restriction of form to the boundary, and the pull back is often omitted, 
for example in Stoke's theorem. In Theorem ll.ll our Dirichlet data Mja^ is the restric- 
tion to the boundary of a form, but in the sense of considering only base points on 
dM and no projection of the fibre on to the tangential component. This Dirichlet data 
together with the integrals of m on a basis of the relative homology group /4(M,3M) 
determines a unique solution to Am = 0. For simplicity we will assume that the said 
integrals are specified to be zero. We note that for the case k — one simply specifies 
the integral of u on each connected component of M with a non-empty boundary. The 
natural Neumann data 1 1 1 is the specification of /* * du and i*5u. Note that as in the 
case k — Q where Neumann data /* * du must have zero integral on the boundary, there 
are compatibility conditions for Neumann data yj . A natural Dirichlet-to-Neumann 
mapping is therefore Ug : D!'[dM) x n"-'^(8M) a''~^-\dM) x Ol'^^dM) given by 

ng(/T,/v) = {i**du,i*8u) 

where Am = 0, i*u = /x, i* *u — /v. Here we use v and T as labels for the normal and 
tangential components in the sense of the domain and range of Ylg. We can now recast 
Theorem ll.ll in terms of this natural data. 

Corollary 2.1. Let M be a smooth compact orientable Riemannian manifold with 
boundary, with dim(M) = m > 2 and metric g, and let k be an integer < k <n. (i) 
The natural Dirichlet-to-Neumann map 11^ /or the k-form Laplace's equation defined 
above is a classical pseudo-differential operator of order L 

( ii) The full symbol ofllg determines the Taylor series ( at the boundary) of the metric in 
boundary normal coordinates. Furthermore for k ^ {(n — 2) /2, (n — l)/2,n} only the 
full symbol of the tangential partTlgt-^ is needed and for and k ^ {0, (« + l)/2, (« + 2)/2} 
only the full symbol of the normal part Tlg^^ is needed. 

The proof of Corollarv l2. ll foUows the proof of Theorem ll.ll in Section|3l 

In the case of electromagnetics note that when formulated in differential forms 
Maxwell's equations ( ll.3> -(1.4) are independent of the ambient Euclidean metric and 
thus invariantly defined. In the inverse boundary value problems for time harmonic 
Maxwell's equations one typically specifies the boundary data invariantly as the tan- 
gential component of E and H. The isotropic case where both the electric and magnetic 
metrics are conformally flat, has been studied by II II and f4l. 

By contrast in elasticity strain is a measure of the distortion of the Euclidean metric, 
and one seeks the elastic tensor with the ambient metric given. This problem is not 
diffeomorphism invariant. 

3 Factorization and symbol calculation 

We consider metrics to be equivalent if they are related by a diffeomorphism which 
fixes points on the boundary. Without loss of generality, therefore, we can assume that 
x„ is the boundary normal coordinate for both metrics. Later we will make a more 
specific choice for the coordinate chart on the boundary. 

We use notation from (4), in particular W)0"''' denotes families of pseudo-differential 
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operators, , in x' such that the /' term in the total symbol vanishes to order I — j at 
x„ = 0, and DO'" '^ is the class of such differential operators. 

LetA'= I . We have that 

A = (Dl +A')I + ED,„ +H{x,D,,), 

where H is a first order system in x' on the bundle of fc-forms and £ is a smooth 
endomorphism of the bundle of A;-forms. We use the notation |^| = y/g{t„t,) for a 
covector ^. 

Proposition 3.1. There exists a B{x,D^,) e TOO'i(3M;R;i2' (M)) such that Oi (B) = 
l^'lJand 

A = iD,J + E + iB)iD,J-iB), 
modulo smoothing and B is unique modulo smoothing. 

Proof. If we expand, we obtain 

dIi + b'-+ED,„ - iEB + i[D„„B]. 

Taking the principal symbol of B as specified we obtain an error, in IDO 'j . Now sup- 
pose we have chosen Bj such that the error, Fj, is in ^HDO'j ^ . Let Bj+i = Bj + C with 
C in W)Oj,[^ . Upon expanding we then obtain an extra term CBj + BjC + Ej with Ej of 
order — /'. Taking a_j(C) — -~j\t,'\^^ai-j{Fj). We have achieved an error one order 
better. Inducting and summing, we achieve an error in MDO^Jj"". 

As the choice at each stage was forced, B is unique. □ 

The importance of this factorization is that B{Q,D^) is equal modulo smooth- 
ing terms to Ag. We will summarise the argument which is identical to that given 
by [6i for the 0-form case. Given a harmonic A:-form u, let v = {Dx„I — iB)u so that 
{Dx„I + E + iB)v = 0. These are both generalised heat equations, the second with 
'time' reversed. As both are smoothing we see that 3„m = Bm + smooth terms and 
hence Ag = B mod W)0^°°. This proves Theorem ll.ll part (i) and part (ii) will follow 
if we can show that two metrics g\,g2 with identical full symbols of B at the boundary 
must agree to infinite order at the boundary. Rather than calculating the full symbol 
of B, we use the calculus of pseudo-differential operators parameterised by the normal 
distance. The advantage is that we need only calculate principal symbols. 

We want to compare the Laplacians associated to two different metrics assumed 
equal up to order I in the normal coordinate. Of course it is immediate from the prin- 
cipal symbol of B that the metrics agree on the boundary so we can take gi — g2— x^k 
for some I > 

First we compare the Hodge star operators. By definition co A *C0 = g{(0, (o)/j and 
we see that 

*i = *2 + x'„a 
where a is a smooth homomorphism from Q.'^ to Q."^'^. 
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Lemma 3.1. If ^^j is the Laplacian on k-forms associated with gj then 

A2-A1 =x'-'FD,^+A 
where F is a smooth endomorphism and A E DO^''. 

Proof. The Laplacian is defined by *d *d and d *d* where d is independent of the 
metric and *2 = *i OC. 

The resuh follows simply by observing that in d*2d*2 and *2d *2 d terms not in Ai 
will vanish to order Z at x„ = unless d is applied to the term. If d is applied once 
to such a term we obtain a first order differential operator vanishing to order I — 1 and 
if twice a zeroth order operator vanishing to order 1 — 2. This is the result desired — 
we know there are no second order terms in Dx^ from our expression for the principal 
symbol. □ 



Lemma 3.2. Let Aj be factored as in Prop \3.1\ with Ej,Bj the corresponding terms. 
We then have that, 

B2-B1 elDO''. 

Proof. Let C = B2 — . As the principal symbols of Bi , B2 agree at jc,, = so we have 
that C is in TOO' '. 

Note that E2 = E\ +x^^'F. Expanding the factorizations for A2, Ai and subtracting 
we have that, 

A2 - Ai === i[C,Dx,]+BiC + CBi +x'-^FDx„ - x'-^F^B + C) +C^. 

After cancelling the x'jj'^FDx,, , we have that 

i[C,Dx„] +B1C + CB1 - ix'-^F{B + C) + de TOO^ '. 

If C e TOO' with 1 <r <l then we have that [C,D,„] e TOO''"', BiC + CBi S 
TOO^'", e TOO^'^", x'-UFB e too' '"', and x'-UFC G TO0'''+'"'. Taking the 
residue modulo TOjO^'*^ we have that i[C,Dx„] +BiC + CBi is congruent to zero modulo 
TOO^'". Recall that 

with Cj G TO)0'[ ^ . The only term of second order is CiZ? + BCi so we deduce that the 
principal symbol of Ci vanishes at x„ = 0. Let cj denote the principal symbol of Cj at 
x„ = 0. We then have that 

(r-;>,- + 2|^'|,c,_i=0 

for each j. Iterating we conclude that cj = for each /' which proves that C G TO)0 ' ^ 
Repeating, the result follows. □ 

It follows from this Lemma that B2—B1 restricted to ;ic„ = is C/ a pseudo-differential 
operator of order 1 — / . Our main result will follow if we can compute the principal 
symbol of this operator and show that it determines k, the lead term of — g2- 



Now 



A2 - Ai = xlP2+x'-'Pi +x'-^Po+x'-'FDx„ 
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with Pj a differential operator in x' of order j and we know from our principal symbol 
computation that P2 is equal to Y,hjDxjDxj where k = —h~^kh~^ where g = dx^ + 

h{x',dx') + 0{x„). 

Arguing as above with C = B2— B\ we have 
j[C,D,„] +B\C + CBi-ix^^^F{B + C)+C'^ =x^^P2+x^n^Pi+x'^^P(i. (3.1) 
So modulo 'M)0^'' we have, 

i[C,Dx,] +B1C + CB1 - ix'-^FBi = x'„P2 +xi-^Pi +x'-W (3.2) 

Let C = Xn ''Cj with Cj e WO j ■' , let cj be the principal symbol of Cj. We have, 

j=o 

i,j<n 

2\^'Uci+lco =Oi(Pi)(^') + /|^'|.F, 
2\i,%C2 + {l-l)ci =ao(Po), 
2g%C2+j + l)ci+; =0, for 1< 7 < / - 2. 

We therefore deduce that 

ij<n 

(3.3) 

where Ki,Li,Mi are computable non-zero constants. (Of course, Mo,Mi = 0. ) 
We want to show that c/ = imphes that kij = for all j, 7. 
If the principal symbol c; = 0, we have taking any component rr of the symbol that 

Kll^Z'^iji^Mj+Llg'lMPlMxA') + gf{iLlFrr + ClCoiPo)rr) = 0. 
ij 

As Pi is a differential operator Gi{Pi)rr{x,t,') is linear in ^' and the final two terms are 
independent of t,'. Since the middle term is not smooth, unless zero, as ^' ^ and the 
other terms are smooth, we deduce that ai (Pi ) = and so we have, 

Kl £ Y.''ij{x)^i^j+\S,f{iLiFrr + CiaoiPo)rr)=0. 

ij<n 

This shows that kij must be a scalar multiple of the identity matrix. Or more invariantly 
that k{x) is a scalar multiple of h{x). To see that k is actually zero, we need to compute 
more precisely. 

For convenience we now reduce to the case of a EucUdean background metric. We 
first prove 

Lemma 3.3. Let 

gi =h+Xnm, 
g2=h + Xnm+x'„r, 

and let *j be the Hodge * operator associated to gj. We then have that *i — *2 modulo 
o(x^+') is independent of m. 
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Proof. Let be a volume form of h and iij a volume form for gj . We have by definition 

1 and g j{ii j,iij) = 1. We then have that/i2 = (1 +x„m(jLi,iLi) +x'„r(jLi,iu))~^^^iLi 
and similarly for . It is now clear that 

1 

V 1 +Xnm{n,n) +xi'r{fi,iu) 

and 



ForanyV,(0ei2*(M) 



^l+Xnm{ju,n) 



(/i+x„OT+4y)(v,«)) {h+x„m){v,()i) 

^yl+x„m{|n,|n) +x'^r(p,ij) ^l+Xnm(ju,fij 

□ 



Lemma 3.4. Let g2 = <ix^ +)^„r + m, and gi = dx^ + m, where m vanishes to second 

order at the origin. Let *j be the Hodge * operator of gj. We then have that *2 — *i 
is independent ofm modulo terms of the form +xj+^w where t vanishes to second 
order at the origin and w is smooth. 

Proof. Let n be the volume form for dx^ . Arguing as above we have that, 

V A (*2 - *i)a) = ^d-' + '-+-nr){yM idx^ + m){y,c,) 
+m{n,n) + x\,r{n,n) ^J\+m{ji,n) 

which modulo terms vanishing appropriately at x = equals the bilinear form 

{dx^+m+x^„r){\ - i(m(/i,/i) +x^„r{p,n)) - {dx^ + m){\ - im(//,/i)), (3.5) 

which upon expanding modulo appropriately vanishing terms equals, 

^n^ -\^nr{n,n)dx^ , (3.6) 

which does not involve m. □ 

Now fix a point p where we will calculate the principal symbol of the difference of 
the Dirichlet to Neumann maps and show that it being zero implies that the next term of 
the difference of the metrics also vanishes there. We take geodesic normal coordinates 
about p in the boundary and then extend normally with respect to . Now we fix a 
point on the boundary p where we will calculate the principal symbol of the difference 
of the Dirichlet to Neumann maps and show that it being zero implies that the next 
term of the difference of the metrics also vanishes there. We choose x! to be a Riemann 



VA(*2-*l)C0 = 

which does not involve m. 
We also have, 
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normal coordinate system on the boundary and then extend normally with respect to 
gi . In particular, we have 

g ~ dx^ + h{x,dx'), 
and on the boundary h{{x\0),dx') = dx'^ + {x'^). 
We then have, using Lemma l3r4l that. 

gi — dx^ + m, 

and 

g2 = dx^ +m+x'„r 

where m = x'j with t vanishing to second order at x' — 0. Let g^ — dx^, and g2 — 
dx^ +x'„r. We then have that *j — *j = x„aj + P; for j = 1,2 with vanishing at ;c' = 
to second order. We also have from our lemmas that 

(*l-*2)-(*l-*2)=-'^«Y 

with y vanishing to second order at x' = 0. 

It is now clear that when computing the lead term of d *2d{*2 — *i) and its appro- 
priate permutations at the point p that we can replace *j by *j, without changing the 

value. So to finish our theorem we take gi — dx^ and g2 ~ dx^ + (1 + x\jX{x'))dx'^ , 
with X a smooth function. 

Let us consider the action of *2 on normal and tangential forms. For a multi 
index / — (/i ,/2, . . . ,4) we use the convention uidxj = Ui^i^^ii^dxi^ Adx,^ A...iix,-^. 
We will denote by /' the complimentary multi index with /' = (yi, • • • ,jn-k) where 
(i'l, . . . ji, . . . ,jn-k) is an even permutation of 1, ... ,n. 

For the metric .... .3,,) is an oriented orthonormal frame and {dxi, . . . ,dxn) 

an oriented orthonormal co-frame. We get an orthonormal frame for g2 by dividing 
each 3,- by (1 + A,(x')x^)'/^ except for d„. Similarly r]/ = (1 + A,(ji:')ji:^)'/^t/x,- for / < n 
and r[n = dxn is an orthonormal coframe for g2- Applying *2 to a typical normal basis 
A;-form dxj (with n £ /), we have 

; k— 1 

*2{dx,) =(1 +X{x')x'„y^ *2 (ri/) 
=(l+Mx'K)-^(il/0 
={l+X{x')x'„)'^-''{dxj,). 

For a typical basis element of i2f (3M), we have n ^ I and 

*2{dxj) ={l +X{x')xl)-^ *2{^I) 

=il+X{x')xi,)'^-\dx,,). 

Modulo xl^^ , we have that on normal ^-forms (*2 — *i) — — ^) A,*i and on 

tangential A;-forms (*2 ^ *i) =x'„ — ^) A, *i . 

To prove the theorem we only have to consider the action of the difference of the 
operators on a particular ^-form ujdxj. 
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We first establish that there is no contribution from 



*2d*2d — *id*id = (*2 — *{)d *2d + *id{*2 — *i)d. 

Any first order term in the first term of the RHS will be in DO^'' and therefore not 
contribute to the principal symbol. For the second term in the RHS, we have, 

*id{*2 — *i)d{uidxi) = *id ^iiM^') *i ^ ^^dxj Adxi^ 

plus terms involving which will not contribute. Now consider the second d, if it 
applies to uj we get a second order term which we already understand; if it applies to 
Xfi we get a term involving dx„ which will then have no dxn component on applying 
*i again. Thus we get no sufficiently low order contribution to either the zeroth order 
term or to the coefficient of Dx„ . 

Moving on to 

d*2d*2 —d *i d*i = d{*2 — *i)d*2 +d*i d{*2 — *i), 

looking at the first term on the RHS we can equally compute with d{*2 — *i)<i*i as the 
difference will be in a non-contributory residue class. 

First considering the normal case n e / we compute 
d*iu= ^ dxj hdxji. 

Now consider 

d(*2 — *i){d jUidxj Adxji) 
'n±\ 



2 

nil 



kj d ^x^X3yM/ *i {dxj Adxp)^ 

kj (^Ix'^^XdjUidxn+xJdCkdjUifj A*i{dxj Adxp) 



Where the + holds for j = n and the — otherwise, but for this to have a dxj component 
we must have j — n. The second term in the final bracket has a coefficient x'„ so can 
only contribute to the second order part which we already understand. We do have a 
contribution to F from 

— k] Ixlj' ' Xd„uidxi . (3.7) 



We are left with the contribution of<i*i<i(*2 — *i).Ifwe apply <ito(*2 — idxj ) , 
we can drop the terms where d falls on a ^ as these are in lOO^''. So on applying 

d{*2 — *\) lo u/dx/ we are left with. 



n + 1 
2 

On applying *i , we get 



xl^X^djUjdxj Adxp + {I — l)xlj ^XuidxnAdxj/ 1 . 
V ye/ J 



^xJ,X^8yM/c/x/^ + - 1)4 ^Xuidxi^ 
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where Ij is simply / with j deleted. 

We now apply d again to get the final contributions. We get another contribution to 
F identical to Eauation l3.7l so that 

and an x^^^m/ which gives 

Substituting back in to Eauation l3.3l we see that if c; = then kij = and Theorem ll.il 
is proved for the case for k {n+ l)/2. Now we consider the case of tangential data 
that is uidx] with n ^ I. 

A similar argument applies to a tangential form, with k ^ (n — 1 )/2 and Theo- 
rem ll.ll is proved. 

Proof of Corollan/ M. 1\ First observe that given the induced metric on the boundary the 
data {i*u, i* * u) determines M|aM- In a neighbourhood's of the boundary a A;-form u can 
be expressed as 

u— ^ ujdxj + ^ U(^f, j'^dx„ /\dxj 

\I\=k,n0 \J\=k-l.n^J 

SO that 

dfi-idu— ^ 3„M/ dxi 

\I\=k,n0 

and 

This shows that K-Jlg is a pseudo differential operator of order 1. Notice now that for a 
harmonic ^-form m, v = *m is a harmonic n — k form for which both the tangential and 
normal parts of Dirichlet and Neumann data are exchanged i*v — i* *u, i* *v — ±i*u, 
i* *dv = ±i*5u and /*5v — ±i*du. It follows that K^Tlg is also a pseudo differential 
operator of order 1. Thus we have proved part (i) of Corollarv l2.1l For part (ii) notice 
that for k = the result is proved by |6l, and here du — identically so the normal 
part of the Neumann data gives us no information. Similarly for k = n the tangential 
part of the Neumann data vanishes. For the case < k < n and k ^ (n ~ l)/2 where 
any tangential-tangential diagonal component of Ag determines the Taylor series we 
require, it is clear that the Taylor series is determined as long as we have We have 
the principal symbol Oi (ngxx)(^) — *d\^\g ^"d to finish the proof for this case we show 
that this determines *a at each point on the boundary. 

Fix a point on dM and choose any multi-indices Iq and Jo such that go{t„t,) :— 
(ai(ngxT)(^)/oio)^ is ^ non-zero quadratic function of ^. Now g{^,^) — ago{^,^) 
where a — l/(*/o/o)^ is to be determined. Let *oa be the Hodge star on fc-forms 
on the boundary determined by go then =1=3 = a'^^("^''/^*oa and (ai (ngxT)(^)/oyo — 
fyk-(n-2)/i l^lgg. As ^0 is known a is determined provided k^ {n~2) /2, hence we 
have g at the boundary and 

For the case A: ^ {0, (n + 1 )/2, (n + 2)/2} we simply apply the above argument to 
boundary data for *m. As one of the conditions on k must hold it is certainly true that 
the full symbol for the complete Ylg determines the Taylor series, at the boundary, of g. 



=iii-i)('^-ky 
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